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Abstract 

The Navier-Stokes equation driven by heat conduction is studied. As a prototype 
we consider Rayleigh-Benard convection, in the Boussinesq approximation. Under 
a large aspect ratio assumption, which is the case in Rayleigh-Benard experiments 
with Prandtl number close to one, we prove the existence of a global strong solution 
to the 3D Navier-Stokes equation coupled with a heat equation, and the existence 
of a maximal B-attractor. A rigorous two-scale limit is obtained by homogenization 
theory. The mean velocity field is obtained by averaging the two-scale limit over the 
unit torus in the local variable. 

1 Introduction 

In this paper we study the Navier-Stokes equation driven by heat conduction. Under a 
large aspect ratio assumption (the spatial domain being a thin layer) we prove the ex- 
istence of a global strong solution. The aspect ratio T is the width/height ratio of the 
experimental apparatus and T > 40 is considered large, see Hu et al. ||l], ||. For existence 
results for the Navier-Stokes equation in thin domains we particularly refer to Raugel 
and the references therein. Our contribution is the addition of a heat equation where 
the heat conduction is driving the fluid. The prototype we have in mind is the classical 
Benard problem and this work is motivated by the instabilities of roll-patterns observed 
in experiments and simulations. Rayleigh-Benard convection is a model for pattern for- 
mation and has been extensively studied, Busse and Clever Q |5| established the stability 
of straight parallel convection rolls and used them to explain many experimental observa- 
tions, for large Prandtl numbers P = ^ (u is the kinematic viscosity and k is the thermal 
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diffusitivity) . For low Prandtl numbers, the situation is much more complicated and it 
has long been recog nized ||, 0, |, i that mean flows are crucial in understanding the 
complex pattern dynamics observed in experiments. In this paper we establish that this 
complexity is not caused by singularity formation but our ultimate goal is a theoretical 
understanding and a quantitative capture of the mean flow. It is known that wave-number 
distortion, roll curvature and the mean flow make straight convection rolls become unsta- 
ble H, 0, m, m. These effects have b een successfully modelled by Decker and Pesch M 
and simulated. However, the resulting equations contain non-local terms due to the mean- 
flow and are theoretically intractable. It is our hope that our results will put the study 
of the contribution of the mean flow on a rigorous mathematical footing, both simplifying 
and reaching a better theoretical understanding in the process. The experimental difficul- 
ties of measuring weak global flow in the presence of dominant local roll circulations are 
formidable, |^, and a theoretical insight may be crucial in the case of dynamical patterns. 

As a preparation, we prove a new statement of the classical, see Leray |14], existence 
of a global strong solution and a global attractor for the three-dimensional Navier-Stokes 
equation under smallness assumptions on the data, see Ladyshenskaya [15|. For different 
results with large forcing compare Foias and Temam [16| and Sell [17|. This new statement 
and new proof of the theorem are crucial in the statement and the proof of the existence 
of a global solutions of the Rayleigh-Benard problem with a large aspect ratio. The 
existence is proven after an initial time-interval, corresponding to a settling-down period 
in experiments. In experiments in gases (CO2) this settling-down time is a few hours for 
experiments that take a few days, 0. 

We will also prove that the Rayleigh-Benard problem has a global attractor. Then 
theorems of Milnor |jl^, Birnir and Grauer [19|, and Birnir 2C] are used to prove that 
the Rayleigh-Benard problem has a unique maximal B-attractor. This attractor has the 



property that every point attracts a set (of functions) that is not shy |21|, or of positive 
infinite-dimensional measure. The discovery of a spiral-defect chaotic attractor 



1% 11 



23|, in a parameter region were previously only straight rolls were known to be stable 
p^ , p^ j, was one of the more startling results in recent pattern formation theory. The 
experimentally observed attractors are B-attractors, and the spiral-defect chaotic attractor 
and the straight roll attractor may be two minimal B-attractors of the unique maximal 
B-attractor whose existence we prove. 



1.1 The Boussinesq equations 

We start with the Boussinesq equations which are two coupled equations for the fluid 
velocity u, the pressure p and the temperature T, 



du 

'm 

dT 
'dt 
divM 



+ {u ■ V)n - i/An + Vp = ga{T - T2) 

+ {u ■ v)r - kAt = 0, 
0, 



X eQ, t e R+. 



(1.1) 



Here g is the gravitational acceleration and a is the volume expansion coefficient of the 
fluid. Moreover i' and k are the viscosity and conductivity coefficients which determine 
the dimensionless Rayleigh R = ^9{Ti-T2)h Prandtl numbers. We assume that Vl is 
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a rectangular box and fix the temperatures at the bottom Ti and top T2 of the box. The 
box is heated from below so Ti > T2. We can impose periodic boundary conditions for u 
on the lateral sides of the box. However, as the temperature Ti increases these have to 
be relaxed due to the presence of a boundary layer. The velocity is assumed to vanish on 
the horizontal surfaces of the box. Finally, we must supply the appropriate initial data. 



1.2 The homogenization 

The mathematical framework for the homogenization starts by the introduction of a 
(small) parameter e > and a scaling of the Navier-Stokes system above 

du 



-gf + iue- V)r, - e^/^AT, = 0, 
, divue = 0, 



(1.2) 



The bulk of the paper will be devoted to proving that there exist unique functions uq, po 
and To such that 



no, Pe Po, To, 



as e — > in the appropriate Sobolev spaces. At a first glance, sending e to zero seems 
to give the Euler equation, but this is wrong. The limit obtained is viscid and the above 
equation only makes sense for e > 0. However, for any fixed e > the solutions of 
the scaled equations have global existence in two dimensions and the equations possess a 
smooth global attractor in dimensions two or three, see Foias et al. Ladyshenskaya 



1 15, 27] and Sell |17]. The equations for the leading order coefficients, in a power series in 
e, are still evaluated at e = and turn out to be the Navier-Stokes system 



— + (uq ■ Vy)uo - uAyUo + Vypi = ga(T2 - Tq) - VxPo, 
dT 

— + (no • Vj^)To - kAj^To = 0, 



(1.3) 



divy no = 0, divx( / mdy) = 0, 



where x &Vt, y £ T" and r G R"^. Here y = x/eis the local spatial variable and r = t/^/e 
is the scaled fast time variable. T", the unit torus in y, is what is referred to as the unit 
cell in the terminology of homogenization. The initial data is so highly oscillatory that we 
can assume that the boundary conditions are periodic in the local variable y. The Navier- 
Stokes system (1.3) differs from the original system (1.2) in that it has an additional 
forcing term — V^po- This is the obvious influence of the global pressure on the local 
flow. The solutions of the system (1.3) enjoy global (in r) existence in two dimensions and 
the equations possess a smooth (in y) global attractor in dimensions two or three. This, 
possibly high-dimensional, attractor of the local flow is the physically relevant quantity 
for the local flow, except in the strong turbulence limit, when long transients may play a 
role. 
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We will show that the solution of (1.3) is the unique two-scales limit of a sequence of 
solutions to the scaled system (1.2). This uses the weak sequential compactness property 
of reflexive Banach spaces and says that any bounded sequence {u^} in say L^(r2) contains 
a subsequence, still denoted by {u^}, such that for smooth test functions ip{x,y), periodic 
in y 

Ue{x)(p{x,-)dx ^ / / uo{x,y)(p{x,y)dydx. 



n JQ Jt 



The main result (Theorem 7.1) is that if {u^} is a sequence of solutions to the Navier- 
Stokes system (1.2), then the so called two-scales limit uq is the solution to the local 
Navier-Stokes system (1.3). Our proof is based upon a compactness result which was first 
proved by Nguetseng [p^] and then further developed by Allaire [29, pH, 31 1. Moreover, if 



Uq is a globally defined unique solution of the system (1.3), which is the case if uq lies on 
the attractor of (1.3), then, by uniqueness, the whole sequence {u^} two-scale converges 
to Uq. 

The mean field turns out to be 



t 



uo{x,-j=)=l (7r(e„6'o)) (x,s)ds, 





where e„, n = 2, 3, is the unit vector in the vertical direction and 7r(e„^o) denotes the 
projection onto the divergence free part of e„0O) 

vr(e„^o) = -V X {A'^V x ejQ)). 

The mean field is derived from the local Navier-Stokes system (1.3). It gives the con- 
tribution of the conduction to the small scale flow. We have averaged (in y), denoted 
by overbar, over the unit cell T", n = 2, 3. Once we have the local problem (4.3) the 
boundary conditions on the local cell can also be relaxed to capture the contribution of 
(global) convection to the mean field. The mean field with the influence of the convection 
taken into account, turns out not surprisingly to satisfy a forced Euler's equation 

— + Uq ■ Vuq + Vpi = 7r(e„^o), 
where r = t/^/e and V = Vy, with y = x/e. 



1.3 Problem setting 

We let be a rectangular box, of thickness h and Lebesgue measure m(J7), in R*^, n = 2 
or 3. By (cj), z = 1, 2 or z = 1, 2, 3, we denote the canonical basis in or R^, respectively. 
The system (1.2) is equipped with the following initial data: 

n(x,0) = uq{x) and T(x,0) =Tq{x), 
which are assumed to belong to L'^{Q), and boundary data: 

u = at Xn = and at Xn = h. 
As above 

T = Ti at Xn = and T = T2 at x„ = h. 
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Moreover we assume that 



u, Vu, T, VT and p 

are periodic with period / in the horizontal xi-direction, in the two-dimensional case and 
periodic with period / in the horizontal xi-direction and period L in the X2-direction in 
the three-dimensional case. In fact we will without loss of generality assume that / = L 
throughout the paper. This determines the pressure p up to a constant that can be fixed 
by normalization, see Remark 2. 



Since we are interested in the fluctuations in the temperature we follow |26] and put 
e = {T-Ti-^{T2-T,)) 



and 



p = p- ga{xn + ^)(^i - ^2). 
We get the following system which is equivalent to ( |1.1D . 

i/An -|- Vp = gaSnO, 



du , 

- + (..V)n 
89 



divM 



(n • V)^ 
= 0, 



h 



X G 0, t e K^. 



(1.4) 



For the temperature 9 we get initial data 9{x,0) = 9q[x) and the new boundary data 
= at x„ = and at x„ = h. The initial and boundary data for u remain unchanged. 
Moreover, u and 9 and their gradients and p are periodic as above. We will find it useful 
to work with the system in the form (|l.l| ) in some situations, whereas the form (|l.4| ) is 
more suitable in other situations. 



2 The Navier-Stokes equation 

We will denote by | • | and || • || the usual norms in L^(J7) and H^{^}). We denote by || • II2 
the i/^(r2)-norm and by | • |oo the L°°(r2)-norm. Further, |ii|2,oo = 6ss sup|ii(t)|, where 
supremum is taken over all i > 0. 

Let us consider the Navier-Stokes equation for incompressible fluids 

' du 

— + (u- V)n - i/An + Vp = f, ^ , , 

dt ' xGl], iGR+, (2.1) 

div n = 0, 

where u is the velocity, p is the pressure and v is the viscosity, with initial condition 

n(x, 0) = no(x) 

and vanishing boundary conditions on dVt (periodic boundary conditions with mean zero, 
if = T^). f denotes the forcing and Ai is the smallest eigenvalue of —A on with 
vanishing boundary conditions on dQ. We start by an estimate which goes back to Leray 



|14|. For the readers convenience we present the (old) proof since the arguments therein 



will be used repeatedly both in the proof of Theorem 2.2 and Theorem 3.3 below. 
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Lemma 2.1. Every weak solution u to the Navier- Stokes equation (2.1) satisfies the es- 
timate 



Moreover, there exists a sequence tj — ^ oo such that 
\Htj)\\ <3^^. 

Proof. We take the inner product of (2.1) with u and integrate over fl. By the divergence 
theorem and the incompressibihty, we are left with 

~Ht)\' + u\Vuit)\' = l^fit).u{t)dx. 

The Schwarz and Poincare inequalities give 

~\u{t)\' + x,uHt)\' <\f{t)Mt)\, 

so by cancellation of \u{t)\, 

^^\u{t)\+X,u\u{t)\<\f{t)\. 
An integration over (0,t), taking sup in t of /, gives 

\u{t)\ < luole-^'"* + - e-^i'^*). (2.2) 

Next we integrate the inequality 

~Ht)\' + u\vu{tr <\f{t)Mt)\ 

over the interval [ti , t2\ , 

u \Vu{t)\' dt < l{\u{t,)\' - \u{t2)\') + \f{t)\\u{t)\dt. 
Jti ^ Jti 

By (2.2) we get 

u r \Vu{t)\' dt < |no|(^ + |/|2,oo)(e-'^^*^ + e-'^'^'^) + + {t^ - t,)). 

Finally, we let t2 = ti + 1, and choose t sufficiently large to get 



ti 



This implies that there is a set of positive Lebesgue measure in every interval [ti ,ti + 1] 
such that 

for t in this set. I 
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We continue by stating a local existence theorem. 
Theorem 2.1. Suppose that the pressure is normalized, i.e., 

pdx = 0, 



then there exists a unique local solution {u, p) of (2.1) in C([0,t]; [H^^ {Q))'^^'^) , n = 2, 3, 
with initial data no = n(x,0) in {H^{Q))'"'. The local existence time t depends only on the 
Lp'-norms of cwiIuq. 

Kreiss and Lorentz |32| can be consulted for details of the proof. 

It is well-known that global solutions of the Navier-Stokes equations exist, 
u G C(R^; in two dimensions. In three dimensions that analogous statement 

is open. However, if transients are allowed to settle for a sufficiently long time, then global 
solutions exist in three dimensions, after this settling of the initial velocity. 

We now state and prove a result saying that every weak global solution to the three- 
dimensional Navier-Stokes equation becomes a strong solution after some finite time to > 0. 
This also goes back to Leray The existence of a global attractor is due to Ladyshen- 



skaya ||l^, 27| and more recently Temam [p3| . Sell |17] and P. L. Lions [^4[. Both the 
statement and the proof of Theorem 2.2 below are new and they are crucial in the state- 
ment and proof of the new existence Theorem 3.3 below, for the Rayleigh-Benard problem 
in the large aspect ratio. 

Theorem 2.2. Consider the three-dimensional Navier-Stokes equation (2.1). For every 
weak solution Uw G C^(R^; to (2.1) there exists ato < oo, such that if the forc- 

^^9 l/b.oo is sufficiently small there exists a unique strong solution u G C([to, oo); {H'^{^))^) 
to (2.1), with data u{tQ) = n^(to)- Moreover, (2.1) possesses a global attractor. If the 
initial data \uo\ is small, then the solutions have global existence, i.e. to = 0. 

Proof. The subscript w indicates that the solutions are only weakly continuos in t. We 
take the inner product of (2.1) with Au and integrate over fi. Using Schwarz's inequality 
we obtain 

^ ^-|Vu(t)|2 + i/|An(t)|2 = [ f{t)Au{t)dx- I {(uit) ■ V)uit)) ■ Au{t) dx 



'^dt JQ JQ 

< |/(t)||An(t)| + |u(t)-Vn(t)||An(t)|. 

Now 

\u{t) ■ Vu{t)\ < |n(t)|oo|Vn(t)| < K\u{t)\^/^\\u{t)\\l^^\Vu{t)\ 
by the Gagliardo-Nirenberg inequalities, where is a constant. Thus, 

^||Vn(t)|2 ^ ^|^^(^)|2 < (1^(^)1 ^ K\n{t)\'/'\\n{t)\\T\Vnm\An{t)l 
An application of the inequality ab < a} j^v + vl? 12, on the right hand side, gives 

^||Vn(t)|2 + ^|An(t)|2<i(|/(t)|+i^|n(t)|V4||^(i)||3/4|Vn(t)|)2 



<i(|/(t)p + i^VWl'/'ll«Wllf |Vn(t)p) 
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by another apphcation of the inequahty above with u = I. From Lemma 2.1 we get, again 
by the same inequahty, 

^|k*)I' + Hv«(i)P < ^(l/WP + KOI')- 

Adding these inequahties, applying the Sobolev inequality and repeating the use of the 
inequality above results in 

lftMt)f + < cMt)f\Ht)f + C2{\mf + \u{t)\% 

where Ci and C2 are constants. Using Poincare's inequality and Lemma 2.1 we find that 
^|lhWf + ^\Ht)f - 2Cmloo + |no|V^^-*)||.(t)f 



smce 



(|/|2,oo + Wole-'^'^f < 2(|/||^ + Kope-^^i'^*). 

Now the point is that if the coefficient |/|| + |iiope~^^i''* is small, or the forcing is 
small and we have waited a sufficiently long time, to let the initial data \uo\'^e~'^^'^^^ decay, 
then the inequality gives us a bound on ||n(t)||. The argument is as follows. We integrate 
the inequality above over [to,t] to get 

\\u{t)f - AC, A|/li,oo + |uope-2^i''*o)e-^(*-^)||u(s)f d5 

J to 

< ^(l/lloo + Kol V2^-*o)(l - e-/^(*-*o)), 

where (5 = Xiv/2. Now assume that ||'u(s)|| assumes its maximum in the interval [to,t] at 
s = t. Then 

AC 



Now put 



vit) = \\u{t)f, a = ^(|/|i,oo + |uope-2^i^*o) 



and 



M = ^(|/li,oo + K|V2^-*°). 

Then the inequality above can be written as 

v{t) - av'^it) < M. 
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Figure 1. The bound on the norm. 



The graph of the function F{v) = v — av'^ is shown in Figure 1. It is concave and attains 
its maximum at Vmax = l/(4a)^/^. By Lemma 2.1 there exists a to > such that 



vito) < \\u{to)r < 



2 / 3|/|2,oo 



and we can choose |/|2,oo so small that v{to) lies between zero and Vmax, see Figure 1. 
Moreover, v{t) can never reach Vmax, because 

F{v{t)) = v{t) - av\t) <M< Ij^^ = F{Vmax). 

It is clear that we can choose to so large that 

M = ^(|/li,oo + K|V2^-*°) 

3 13 1 

< - 



4 (4a)V3 4 {^\f\l^ + |noPe-2Ai-to)i/3- 

Moreover, notice that the derivative of F in Figure 1 is positive at the point, where F first 
reaches M, i.e., 

F'(v) = 1 - Aav^ = 1 + - 4 = 3>0, 

V V 

which gives the bound 

vit) <— + S 

where 6 is arbitrarily small. The only question left is whether the initial data makes sense 
as a function in H^{Q,). But we have already used above that, by Lemma 2.1, there exists 
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a sequence tj — > oo such that < oo. We now choose the initial time to be the 

smallest tj > to- Then we let this tj be the new tQ. 

Now we apply the local existence Theorem 2.1 and the above bound to get global 
existence. We recall the definition of an absorbing set from Coddington and Levinson 
|35|. A set D C is an absorbing set if for every bounded set M there exists a time t{M) 



such that t > t{M) implies that u{t) G L>, if m(0) E M. The estimate in Lemma 2.1 shows 
that the weak-flow of the weak solution of the Navier-Stokes equations has an absorbing 
set in L^. But we have shown in addition that there exists a ti{M) such that the flow 
is a continuos flow of strong solutions for t > ti{M) and has the absorbing set D C H^. 
Moreover, since H^{Q) is compactly embedded in L^(J7) it follows that the Navier-Stokes 
equation has a global attractor in L'^{Q). In fact one can now show, see [32|, that the 
solutions are spatially smooth and thus D C C°°. It follows, see Hale [^|, Babin and 



Vishik 1 37], Temam |33| and Birnir and Grauer |38], that the Navier-Stokes equation has a 
global attractor consisting of spatially smooth solutions. It is also clear that if |uo| is small 
then v{t) satisfies the above bound and the solutions have global existence for to = 0. ■ 

Remark 1. The statement of Theorem 2.2 is really just a restatement of Leray's classical 
result of global existence for small initial data. The duality in the smallness condition is 
that one can either take small initial data or, with arbitrary initial data, just wait for a 
sufficiently long time. This statement is what most physicist and engineers are interested 
in, at least for low Reynolds numbers. Namely, transients die out and the flow settles 
down to the flow on the attractor in at most a few hours, in experiments, see for instance 



Remark 2. The existence of the pressure term follows from a standard orthogonality 
argument which we present below. Consider the Navier-Stokes equation (^]^) and the 
solution given by Theorem 2.2. We take the inner product by u, integrate over apply 
the divergence theorem, use the incompressibility and collect all terms on one side. This 
gives 



{f — u ■ Vii + vAu) ■ udx = 0. 

Jn cli 



Now we recall that the orthogonal to the divergence free elements in are gradients in 
L^. Therefore, there exists a unique gradient Vp in L^(r2) given by 

du 

V» = / — u ■ Vn -|- z^An. 

ot 

We say that a set M is invariant under the flow, deflned by a nonlinear semi-group 
S{t), if S{t)M = M, i.e. M is both positively and negatively invariant. An attractor A 
is an invariant set which attracts a neighbourhood U or for all £ U, S{t)xo converges 
to „4 as t ^ oo. The largest such set U is called the basin of attraction of the attractor 
A. If the basin of attraction of an attractor A contains all bounded sets of X and A is 
compact, then we say that A is the global attractor. 

We have shown above that the Navier-Stokes equations have a global attractor. But we 
have not explained what the semi- group is that shrinks bounded sets onto the attractor. 



It is possible to define a semi-group on the space of weak-solutions, see Sell |17], but in 
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light of Theorem 2.2 we can do much better. Namely, a bounded subset of a Banach 
space is a complete metric space and we have shown that for every bounded set M C 
eventually the solution starting in M will lie in an absorbing set D C C°°. It is the 
semi-group defined on this complete metric space D C that has the global attractor 
A = uj{D), where uj denotes the w-limit set, see [38|. A is invariant, so on it we can solve 



the Navier-Stokes equation in backward time, it is non-empty and compact and attracts 
a neighbourhood of itself, see |38|, moreover, it also attracts every bounded set in L^. 



This attractor consists of spatially smooth solutions and has finite Hausdorff and fractal 



dimensions, see 40, 41, |4^. The dimension estimates are however large, see |^3[, and 
not necessarily a good indication of the size of the attractor. 

We focus on the core of the attractor A which is called the basic attractor B. For now, 
let us assume that the Banach space X is finite dimensional. An attractor B is called a 
basic attractor if 

• The basin of attraction of B has positive measure. 

• There exists no strictly smaller B' C B, such that up to sets of measure zero, hasin[B) 
Cbasin{B'). 

A global basic attractor is a basic attractor whose basin of attraction contains all bounded 
sets of X up to sets of measure zero. 



A theorem by Milnor |15] states that in finite dimensions there exists a unique decom- 



position of the attractor A, 
A = BUC 

where ;S is a basic attractor and C is a remainder such that m{basin{C)\hasin{B)) = 
where m is the Lebesque measure. 

The notion of a basic attractor can be extended to infinite dimensions and Milnor's 



Theorem was first proven in an infinite-dimensional setting by Birnir and Grauer |19]. 
They used cumbersome projections to a finite dimensional space where the ^-attractor 
resides. A more elegant notion of a ^S-attractor requires an extension of the concepts of 
measure zero and almost everywhere. Their counterparts in infinite dimensions are shy 



and prevalent sets respectively, see Hunt et al. |21|. These are defined in the following 
way. Let X denote a Banach space. We denote hy S + v the translate of the set S <Z X hy 
a vector v G X. A measure ^ is said to be transverse to a Borel set S" C X if the following 
two conditions hold: 

• There exists a compact set U C X for which < fJ,{U) < oo. 

• fi{S + v) = for every v £ X. 

A Borel set S C X is called shy if there exists a compactly supported measure transverse 
to S. More generally, a subset of X is called shy if it is contained in a shy Borel set. The 
complement of a shy set is called a prevalent set. 

The infinite-dimensional analogue of Milnor's Theorem can now be stated. 

Theorem 2.3. Let A be the compact attractor of a continuous map S{t) on a separable 
Banach space X . Then A can be decomposed into a maximal basic attractor B and a 
remainder C, 

A = B\JC 
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such that the realm of attraction ofB is prevalent but the realm of attraction ofC, excluding 
points that are attracted to B, is shy. 



For a proof of this Theorem see Birnir [20| and |43]. It imphes that the Navier-Stokes 



equation has a maximal i3-attractor in two and three dimensions. If B can be decomposed 
into finitely many (disjoint) minimal i3-attr actors, then the union of the realms of attrac- 
tion of these minimal ;B-attractors is the whole space c R))",n = 2,3. The realm 
of attraction is a slight generalization of the basin of attraction, as the basin is an open 
set, but the realm can be either open or closed. 

For systems that are simple enough, the basic attractor contains only the stable tra- 
jectories of the global attractor. When transients are ignored, one will only see the basic 
attractor in physical experiments and numerical simulations, not the remainder C = J^B 
of the attractor. All the relevant dynamics are therefore contained in the basic attractor. 
Ladyshenskaya |4^ gives more examples of B-attractors for the Navier-Stokes equation 
with nonlinear viscosity. 

3 The Rayleigh-Benard convection 

3.1 Existence results 

Now consider the Boussinesq system (1.1) equipped with initial data 

u(x,0) = uq{x) and T(x,0) =Tq{x). 

The boundary conditions are periodic on the vertical surfaces = 0, / for k < n. Fur- 
thermore, on the horizontal surfaces Xn = and Xn = h, the velocity u vanishes (no-slip 
condition) and the temperature T is kept constant, i.e., 

T = Ti on x„ = 0, T = T2 on Xn = h. 

We want to state a global existence theorem for the solution to the system (1.1). For 



that purpose we start with the following maximum principle, c.f. |26]: 
Lemma 3.1. Suppose that u and T solve (1-1) ■ If 

T2 <r(x,o) <ri, 

for a.e. x € il, then 

T2 < T{x,t) < Ti, 
for a.e. x E and all t >0. 
Proof. We consider 

(T — Ti)_(_(x, t) = ess sup(T — Ti)(x, t). 

The second equation in (3.1) gives 

^(T - ri)+ + (n • V)(T - ri)+ - kA{T - T,)+ = 0. 
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We multiply this equation by (T — ri)_|_, integrate over and use the divergence theorem, 
to get 

~\{T- ri)+(t)p + k\v{t - ri)+(i)|2 = 0. 

Thus, by Poincare's inequality 

~\{T- ri)+(i)|' + \in\{T - ri)+(t)p < 0, 

where again Ai is the first eigenvalue of the negative Laplacian with vanishing boundary 
conditions on J7. Consequently 

|(r-ri)+(t)|<|(T-ri)H-(o)|e-^i«*, 

which shows that 

(r-ri)+(.,t) = o, 

for ah t > 0, if 

(r-Ti)+(.,0)=0. 

Similarly wc get that 

{T-T2)-{x,t) =ess sup(-(T-r2)(a:;,t)) = 0, 
sen 

if (T - r2)_(x,0) = and we conclude that 

T2<T{x,t)<T^. ■ 

Remark 3. Lemma 3.1 actually yields a uniform bound on T in L^([0, s] x for any 
s > 0. This also gives a bound for 0, see below. 

Lemma 3.2. Every weak solution u to the Boussinesq equations (1-4) satisfies the esti- 
mate 

Ht)\ < l^ole-^^^* + ^ (1 - e-^i^*), 

mt)\ < — , 

where K = ga{Ti-T2)L/3^/^ . The equations possess an absorbing set in (L^(r2))^, defined 
by 



Ht)\ + \9{t)\ <(— + -r^) Kh^''' + S, 
\ga \\v ) 



where b is arbitrarily small. Moreover, there exists a sequence tj oo such that 

Mtj)f + \\e{t,)f<K,, 

where 

(h ^ {T^-T2? 
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Proof. We recall the relationship between T and 

T = e + Ti-^{Ti-T2). 
The maximum principle in Lemma 3.1, for T, T2 < T < Ti, implies that 

_(1 _ ^)(Ti _ T2) < < ^(Ti - T2). 

Thus 

1^12 < (^1 - T^fL'h/Z. 

Now consider the system (1.4). We multiply the first equation by u and integrate over 
Integration by parts and Schwarz's inequality give 

^|Kt)|2 + Aiz.Ki)|2<H^(t)||«„(t)|. 

Thus by the use of the bound for 9 above, and Poincare's inequality, we get 

j^\u{t)\ + \iu\u{t)\ < ga\0{t)\ < ga{Tr - T2)L{h/3)'/'' . 
Integration in t then gives 

\u{t)\ < luole-^i'^* + ^^(1 - e-^i^*)> 

where K = ga{Ti - r2)V2L/(3)i/2. 

We combine the bounds for 6 and u to get the absorbing set 



1^1 



\ga \xv ) 



in where 

^=g^ 31/2 ^ 

and b is arbitrarily small. 

The last statement of the lemma is proven by a straightforward application of Lemma 
2.1, to the equations (1.4), if we recall that the nonlinear term did not play a role in the 
proof of Lemma 2.1. Namely, for the first equation in (1.4), 

l/lioo=5V|0|2 = K2^, 



and for the second equation 

2 _ (Ti - T^f , ,2_ (Ti-T2) ^K^ 

|/|2,oo- p M - ^ 



We divide these by the decay coefficients A^j^^ and k^, respectively and add them. This 
produces the bound on the norm for the sequence tj. ■ 
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Theorem 3.1. Suppose that the pressure is normalized, i.e., 
pdx = 0, 



m{n) 



n 



then there exists a unique local solution {u,6,p) of (1.4) in C([0,t]; {H^^ {^))'^^'^) , n = 2, 3, 
with initial data {uo,9o){x) in {H^ {Q))"'~^^ . 

Proof. The maximum principle in Lemma 3.1, for T, T2 < T < Ti, imphes a maximum 
principle for 9, 

-(l-^)(Ti-T2)<0<^(Ti-r2), 

by the relationship between T and 9 

T = 9 + Ti- ^{Ti-T2). 
h 

This implies that 

\9\l < (Ti - T2fL'^h/3. 

The rest of the proof, using this bound on 0, is similar to the proof of the local existence 
of the Navier Stokes equation. Kreiss and Lorentz [^] can be consulted for details. Then 
given the local solution {u,9){x,t) the pressure is recovered as in Remark 2. ■ 

In two dimensions we have the following global existence result. 

Theorem 3.2. The Boussinesq system (1-4) has a unique global solution (n, 9) in 
C(R+; (-ff^(r2))^), where Q is a bounded open set in R^. Moreover, the system (I.4) 
possesses a global attractor in 



Proof. A proof of Theorem 3.2 can be found in Foias et. al. |26|. 

In three dimensions the following global existence result holds true: 

Theorem 3.3. Consider the Boussinesq system (1.4)- For every weak solution u^, 9^ 
in Cw(R^',{L'^{^))^), ^ a, bounded open set in with a large aspect ratio, ^ >> 

^'^^uK^'' ^ ' ^^^re L is the width (radius) and h is the height ofCl; there exist a time 
to, such that there exists a unique strong solution (n, 9) in C{[to,oo);(H^{il,))^), t > to, 
with initial data u{to) = 'U^(to) CLi^-d 9{to) = ^«,(to)- Moreover, the system (I.4) possesses 
a global attractor in 

Proof. The subscript w indicates that the weak solutions are only weakly continuos in t. 
We multiply the second equation in (1.4) above by 9 and treat it in the same way as u in 
the proof of Lemma 2.1 to get 

+ ^|Vg(t)p < - \9{t)\\Un{t)l 



and by adding and subtracting f3\9\ and applying Poincare's inequality, we get 

a!' 



^^-if»(*)i + mt)\ + (a;"" - ^)iv9(()i < ' \ " \u„(t)\ 
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< ^^^:^^^ K|e-"^^* + ^(1 - e-'^"^')) < ci + cse-^i"^*. 
Then integrating with respect to t, we get 

„ (p-l3t _ p-\\Kt\ 

Next we multiply the u equation in (1.4) by Au and integrate over f2. By integration 
by parts and Schwarz's inequality 

^^|Vt/(t)|2 + v\^u{t)f < ga\Ve{t)\\Vunit)\ + \u{t) ■ Vu{t)\\Auit)\ 

< 5«|v^(t)||Vt.„(i)| + c(i + Ar' + A^'f/«|«(^)|V4|v«(^)||An(^)^/^ 

by the Gagliardo-Nirenberg inequalities, where we have used that the Sobolev norm 
is bounded by 

||n||2 < C{\u{t)\'' + |Vt^(i)|' + |An(t)p)V2 < c(l + A^' + Ar')'/'|Au(t)|, 
by Poincare's inequality. We use Young's inequality to eliminate lA^I^/^, namely 

so 

~\Vu{t)f + ^|A^f)|2 < gamt)\\Vu{t)\ + Co\u{t)f\Vu{t)f 
and by Poincare's inequality 

||Vn(t)| + ^\Vu{t)\ < ga\Ve{t)\ + Co\u{t)\'\Vu{t)\' . 
Then we integrate the equation 

^|Vn(t)| + l3\Vu{t)\ - CoHt)\''\Vu{t)f < ga\Ve{t)\, 
where (3 = \\u/2. We integrate from the initial time Iq, to get 



\Vu{t)\-Co f \u{s)\^\Vu{s)\'e 
Jto 



<ci + c2e-^(*-*°) + C3e-^^'^(*-*»), 

by the above inequality for J \ V9\dt. Now if |Vn(s)| attains its maximum on the interval 
[to, at s = t, then 



Jto Jto 
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< (|uop(e-3/^*o - e-3/^*)e-/^* + 3^(1 - e-/^(*-*o))) , 

where we have used the bound 

\u{t)\ < |uo|e-^^"* + ^^(1 - e-^i"*), 

from Lemma 3.2 and the inequality ab < a^/2 + 6^/2. Thus, if we put |Vu(t)| = v{t), we 
obtain the inequality 

v{t) - av^it) < M. 

The constants a and M are, 



and 



M = 2ci = 2^\j^ = 2ga{T, - T-,fL/?>y^\,uh"\ 



This means that for a large aspect ratio L/h >> ^" '^'^^uk^^ ^ ' ^ becomes small. Now 
we repeat the arguments from the proof of Theorem 2.2 and conclude that the function 
F{v) = V — av^ is concave and attains its maximum at Vmax = l/(7a)^/^. This maximum 
is 

6 1 



7 (7a)V6 ' 

so that v{t) can not escape beyond its maximum, see Figure 1. Moreover arguing as in 
Theorem 2.2 we conclude that the derivative of F is positive at the point where F first 
reaches M and this gives us the bound 

7M 

"^^^ 

The last step is to get a bound for \V9\. We multiply the 6 equation in (1.4) by A0, to 

get 

1 |jV0(t)p + K|A0(i)P < ^^'~J^\ ve{t)\\Vu{t)\ + \u{t) . V9{t)\\A9{t)\ 
by Schwarz's inequality 

- + \u\e\V9it)k\A9it)\, 
by Holder's inequality 

- ^^^^^|V^(i)l|Vn(t)| + C{\Vuit)\''\9\/S' + S\A9\)\A9{t)\, 

by Poincare's and Sobolev's inequalities, because 
\u\(i < K\\u\\ < C\Vu\, 



Existence and Homogenization of the Rayleigh-Benard Problem 



153 



as above, and because, 

\VO\3 < C\\Ve\\y2 < K{^^ + 5\M\) < + 5\A9\) 
by two applications of interpolation, where 6 is small. Thus 

by Young's inequality. Now moving the |A6'(t)p term over to the left hand side of the 
inequality and applying Poincare's inequality, we get 

~mt)\' + Ai(« - 35)\V9{t)\' < ~mt)\' + - 36)\Ae{t)\' < C3, 

since 

^^'~J'\ ve{t)\\Vn{t)\ < X,6\Ve{t)\' + i^^^I^|Vn(t)p, 

by the inequality ab < + 5^/4. Thus 

\Ve{t)\^ < |V0(to)|'e-T(*-*°) + (2c3/7)(e-^*» - e"^*), 

where 7 = 2{k — 35). 

We can now put the estimates for \Vu\ and 1X76*1 together to get an absorbing set 

\Vu{t)\ + |V6'(t)| < constant + e, 

where e is arbitrarily small for t large enough. Namely, by Lemma 3.2, there exists a tj, 
such that {u,9){tj) £ H^(Q) and the global bound on the norm holds for t = to = tj. 
Combined with the local existence Theorem 3.1, the a priori bound above now gives the 
existence of a global solution and an absorbing set in H^{Q), for t > tQ. The existence of 
a global attractor in then follows from the compact embedding of {H^{Q))'^ in 

see Hale p6|, Babin and Vishik [R^], Temam [pSl and Birnir and Grauer 



Corollary 3.1. The Boussinesq system (I.4) has a maximal B- attractor whose basin is a 
prevalent set in (L2(R")"+1), n = 2, 3. 

The proof is a straight-forward application of Theorem 2.3 

Theorem 3.3 says that in experiments in pattern formation where one has a large aspect 
ratio, it is only necessary to wait a short time to have global spatially smooth solutions. 
The smoothness follows from the smoothness of the nonlinear semigroup, see Kreiss and 



Lorentz [32|. Corollary 3.1 says that what is observed in the experiments after the initial 
settling-down time, is a ;S-attractor that is a component of the maximal ;S-attractor. In 
other word this attractor is observed for an open set of initial conditions. (However, not 
necessarily an open set of parameter values.) It is then useful to know how long one has 
to wait and this time is easily estimated. It is the time it takes the exponentially decaying 
initial data to become comparable in size with the absorbing set. Namely, 



t = -r^ln 



Aii^|no| 



3I/2A1 



ga{T, - T2)Lhy^ 



In experiments for Prandtl numbers close to one, this time is measured in hours for 
experiments that take days S. 
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4 Scaling and expansions 

As a starting point we introduce a small "scaling" parameter e > and consider the scaled 
system 

^ + K • v)r, - 6^«Ar, = 0, X G ^^, t G R+, (4.i) 

at 

, div = 0, 

equipped with initial data 

'Ue(a;, 0) = a{x) and re(a;, 0) = b{x). 

The boundary conditions are periodic on the vertical surfaces Xk = Q,l for k < n. Fur- 
thermore, on the horizontal surfaces x„ = and .x,„ = h, the velocity vanishes (no-slip 
condition) and the temperature is kept constant, i.e., 

= Ti on Xn = 0, = Ta on x„ = h. 

A small value of e corresponds to a high Reynolds number and the solutions will likely 
become turbulent. We also consider the scaled system corresponding to (1.4), i.e., 

( dUe 

f + K.V)«.-.^KA9. = 2i^(„.)„, ««.*<^R^- (4.2) 

divug = 0, 

For the temperature 0^ we get initial data 0e(x,O) = 9'^{x) and the new boundary data 
= at x„ = and at x„ = h. The initial and boundary data for remain unchanged. 

Moreover, and 9^ and their gradients and are periodic as above. 

It turns out that the value 7 = 3/2 is critical. To see that we perform a multiple scales 

expansion technique of the unknown quantities Ug, p^ and 6^ and assume that 



X t 

Pe{x,t) = ^eVi(a;, —,t, -^), 



i=0 



00 

X t 



i=0 

where Ui, pi and 6i are all assumed to be r"-periodic with respect to j/ G R", n = 2, 3, 
T'* being the usual unit torus in R'*. If we put y = x/e^ and r = t/e^, the chain rule 
transforms the differential operators as 

d d I d d did 
dt dt Ot ' dx dx e/* dy ' 



Existence and Homogenization of the Rayleigh-Benard Problem 



155 



The question is: Can we preserve the structure from the original problem in the leading 
order approximation? One sees that if the viscosity and conductivity terms scales like 
e^/^, then the choices fi = 1 and p = (3 = 1/2 preserve all quantities. By substituting 
all this into the system (4.2) we can equate the powers of e. This formal manipulation 
shows that the functions Ui and 6i, i = 0, 1,... are all independent of t and that po is 
independent of y. Formally this means that there are two scales in space and in order to 
preserve the evolutionary behaviour, in the leading order approximation, a chance of time 
scale t I— > ^jy2 = becomes necessary, c.f. Lions [^]. It also follows, from the equations, 
that the functions pi are independent of t. The leading order system reads 

— + (uq ■ Vy)uo - uAyUo + Vypi = gacnOo - VxPo, 



do 

-^ + {UQ- Vy)eO - KAyOo = 0, 

divyuo = 0, divx( / uody) = 0, 



(4.3) 



on $7 X T" X R+ with T"-periodicity in y. 

We see that the system (4.3) differs from the system (1.4) only by the additional forcing 
term coming from the "global" pressure gradient. Thus the local problem (on the e-scale) 
has two pressure gradients, one local VyPi and one global VxPo- 

Remark 4. The scaling discussed above is isotropic but the convection rolls in Rayleigh- 
Benard convection have a preferred direction, say along the X2-axis. This commands the 
use of anisotropic scaling for the amplitude equations and it turns out that the scaling by 
e^/'^ is the one pertaining to the other directions, perpendicular to the orientation of the 
rolls. The details of this affine scaling will be spelled out elsewhere. 



4.1 Existence results 

In three dimensions the following global existence result holds true: 

Theorem 4.1. Consider the Boussinesq system (4-2). For every fixed value of e > the 
following holds true: For every weak solution Uw,e, dw,e in C^(R"''; {LP')^{VL)), VL a hounded 

open set in R^ with a large aspect ratio, ji » ^° '^'^^vk'^^ ^ ' "^here L is the length and h is 
the height ofQ; there exists a time to, such that there exists a unique strong solution u^, 9^ 
in C{[to,oo)]{H'^Y{^l)) , t > to, with initial data Ue{to) = iit«,e(to) O'^d 0^{to) = Ow,e{to)- 
Moreover, the system (4-^) possesses a global attractor in (L^)^($7). 

Proof. Theorem 4.1 follows from Theorem 3.3 if we set the viscosity and heat conductivity 
in Theorem 3.3 equal to e^^^z^ and e'^/^K, respectively. 



Remark 5. The existence of a pressure Pe G C([to, oo); {H^)^^.)/!^) follows by the same 
arguments as those in Remark 2 after the proof of Theorem 2.2. 



We conclude with the existence result for the two-scale system (4.3). 
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Theorem 4.2. Consider the system (4-3) with initial data UQ{x,y,0) = a{x,y) and 
9Q{x,y,0) = b{x,y), with zero mean over T". For almost every x £ Q the system 
has a unique global strong solution uq G C([to, oo); (//""^(T"))"), n = 2,3, and 9q G 
C([to,oo);i?Hr'^)). Moreover, pi G C([to, oo); i7i(r")/R). Integrated over {uo)i, 
9o and pi all belong to L?'{Q). Finally, po G C{[to,oo);H^{^l.)). In two dimensions to = 
and in three dimensions to = to{\a\, \b\) > 0, in general. Moreover for almost every x € 
the system (4-3) possesses a global attractor in (L^(T"))" x lP'[T'^), n = 2, 3. 

The proof in the two-dimensional case is relatively straightforward and details can be 
found in Foias et al. p^]. 

Proof. In the three-dimensional case the proof is similar to the proof of Theorem 3.3 but 
simpler. First we multiply the first equation in (4.3) by 9 and integrate by parts to get 

\9o{t)\ < |6|e-^* 
after integration in t, where 7 = Aik. In 

gae39o - V^Po = gaTT2{e39o), 
where vr2(e3^o) denotes the projection onto the divergence free part of es^o- Hence, 

\gae3e0 - VxPol < 5a|^o| < 5a|^|e"^*. 

This means that for t sufficiently large the forcing in the first equation of (4.3) becomes 
small. Thus we obtain the existence of unique global solution uq G C{[to, 00); (i7^(T"))^) 
to the first equation of (4.3), by Theorem 2.2, and an absorbing set in this space. Since 
uo{t) G {H^{T^))^ we also obtain global existence for in the second equation of (4.3), 
just as in the proof of Theorem 3.3 and the existence of an absorbing set in this space as 
well. The existence of the attractors follows from the fact that the equations possess an 
absorbing set in (ii"i(r"))3 x H^{T'') which is compactly embedded in (^^(r"))^ x L2(r"). 



46 1 we show that 



Remark 6. In Section 7, we prove, by using homogenization theory, that the solutions 
of the scaled system (4.2) converges in the two-scale sense, see [p8| , to the unique solution 
to the system (4.3). 

Remark 7. By averaging the system (4.3) over the unit torus in local variable y we obtain 
the mean-field corresponding to the scaled system (4.2) in Section 8. The understanding 
of the effects of this field on the system is crucial in the theoretical understanding of the 
complex patterns in Rayleigh-Benard convection. 

5 A priori estimates 

We are interested in the asymptotic behaviour of the system (4.2) as e ^ 0. In order to 
accomplish this we need to establish uniform (in e) bounds on u^, Pe and 9^. 

This is problematic, since the forcing term in the Navier-Stokes equation in (4.2) in- 
volves 9e which can be highly oscillatory for small values of e. We begin with the following: 
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Lemma 5.1. Let u, p and 9 be the solution to the system 

^ + (n • V)u — vAu + Vp = CnO, 
ot 

90 An i\ xGfi, iGR+. (5.1) 

+ (u • V)9 - kAO^ = {u)n, ' ^ ^ 

, divu = 0, 

Suppose that u, p and 9 are all T" periodic. If 
[ Ui{y,0)dy= [ 9iy,0)dy = 0, 

then 

I Ui{y,t)dy= 9{y,t)dy = 0, 
for all t>0. 

Proof. By the equation (5.1) we have 
f du f 

/ {-^)idy= / {{en9)i-{{u-V)u)i + iy{Au)i-{Vp)i)dy 



and 



/ (S)^y= / {{u)n-{u-V)9 + nA9)dy. 



All terms involving spatial derivatives will vanish by the r"-periodicity on the unit torus 
and by the incompressibility of u. Moreover, (e„0)i = 0, for i < n, so the system reduces 
to 



di 



/ iu)ndy = / 9dy 



and 



— / 9dy = I (u)ndy. 
dt Jt^ Jt" 



dt 

We solve this system of ODE's with the initial condition 

/ Ui{y,{))dy= [ 9{y,0)dy = 0, 
to get the trivial solution, 

/ Ui{y,t)dy=i 9{y,t)dy = 0, 

Jrpn Jj-n 

for all t > 0. 
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We are clearly making a strong assumption assuming that the solution is periodic in 
the above Lemma. In the applications we have in mind it will not be exactly periodic but 
it will be sufficiently oscillatory so that periodicity on a small scale is a reasonable model. 
Flow in porous media is obviously one example of such a situation but others include, a 
mixture of hot and cold fluid, a two fluid mixture (oil and water), a mixture of snow and 
air, or water and sediment, and a turbulent fluid. This hypothesis is similar to the model 
one uses in homogenization of materials and the test of the hypothesis is how well one 
captures averaged quantities. 

The next lemma is a Poincare inequality, c.f. L. Tartar (Lemma 1 in the Appendix of 
p7|). Our version of the lemma differs from Tartar's in that we do not have a vanishing 
boundary condition on the local tori. Instead we benefit from the result of Lemma 5.1, 
that the mean value vanishes. 

Lemma 5.2. Let ui E -^^^([0, s]; H^{Vl)) and assume that u is periodic and that initially u 
has mean value zero over the unit torus T", then 

u{x,t)\^dxdt<e^C I [ \Vu{x,t)\'^dxdt, (5.2) 
Jn Jo Jn 

where C is a constant independent of e. 

Proof. We apply the Poincare inequality on T" which, by the result of Lemma 5.1, gives 

u{y,t)\'^dy <Ci I \Vyu{y,t)\'^dy. 
A change of variables x = ey yields 



\u{x,t)\^dx < e-^Ci \V:,;u{x,t)\'^dx. 

eT" JtT" 

We note that the constant Ci will be the same for all eT'^-cubes in the interior of 17. For 
eT'^-cubes intersecting 90, u will vanish at, at least, one point and the usual Poincare 
inequality applies. Let C2 denote the maximum of all the constants for these cubes. A 
summation over all eT"-cubes gives 



u{x,t)ydx <Ce^ I |V^M(x,t)rdx, 
n Jn 

where C = maxjCi, C2}. Finally, an integration with respect to t gives the desired 
inequality. ■ 

Lemma 5.3. Let and 6^ satisfy (2.1) and assume that the initial data in (4-2) is 
bounded independent of e, then 

e~"^'^^lke||L2([io,s];L2{C)) < C (5.3) 

and 

e^^^\\Vu,\\L2^[to,s];LHn)) ^ (5-4) 
for any s > 0, where C is independent of e. 
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In the two-dimensional case to = and in the three-dimensional case to > in general. 
This will hold true throughout the paper. 

Proof. Let us consider the first equation in (4.2). We multiply by and integrate over 
f2x]to,s[. By the incompressibility we get, by using the Schwarz's inequality 



< 



1, 



lL2([to,s];L2(n))ll«e||L2([to,s];L2{f7)) +T^l^e(io)| • 



By Lemma 3.1, see also the proof of Theorem 3.2, ||^e||L2([to,s];L2(n)) is bounded, and 
assuming that to is a time where |Vtt| is finite, see Lemma 2.1, we can absorb the term 
2|'Ug(to)P into the time integral using Lemma 5.2. This gives the estimate. 



-\Ue{s)\'^ + e^/'^\\VUe\\l2^^to,s];L2m ^ ^ M L^[to ,s];L^n)) 



2' 

so that, 

e^^^l|VUe||i2([tg,,].i2(n)) < C\\Ue\\Lmto,s];L^{n))- 

Now we recall that G L^{[to, s];H^{n)) by Lemma 3.2 (or Lemma 2.1). Thus, by Lemma 
5.2, 

e^^'^\\VUe\\L2{[to,s];L^{Q)) < C. 

An application of Lemma 5.2 once again gives the desired result, i.e., 

^'^^^\\Ue\\L^([to,s];L^in)) < C 

If the initial data Ue{to) and Oe{to) is bounded, in L^(ri), independent of e then C will also 
be independent of e. ■ 

We continue with a few consequences of the above results: 
Corollary 5.1. Consider the first equation in (4-2). The convection term is hounded, 

||(tXe • V)^ie||L2([to,s];Ll(n)) < (5.5) 

where C is independent of e. 

Proof. By the Schwarz's inequality. Lemma 3.2 and Lemma 5.3 it immediately follows 
that 

\\ue ■ V'Ue||i2([jp^s].^i(n)) < e~-^/^||ug||^2([iQ^s].^2(n))e^/^||Vng||^2([iQ^5].i2(n)) < C^. | 

Corollary 5.2. Consider the first equation in (4-2). The time derivative is hounded, 

ll^llL2([to,s];L2(n)) < C", (5.6) 
and the pressure is hounded, 

\\Pt\\L-^{[to,s\;m{Q.)/-R) < C, (5.7) 

where C is independent of e. 

Proof. (5.6) follows from duality by Corollary 5.1 and (5.7) follows by the Remark 2. ■ 



160 



B. Birnir and N. Svanstedt 



6 Two-scale convergence 

In this section we recaU the technically useful concept of two-scale convergence ([^ and 
Ipsll ), for the case when the functions also depend on a time variable. Let us consider the 
space C°°(T") of smooth periodic functions, with unit period, in R"-. 

Definition 6.1. A sequence {u^} in L'^{[0, s]; is said to two-scale converge to uq = 

uo{x,y,T) in L2([o,s];L2(Q x T")) if, for any ip G (7^(0 x [0, s]; C°°(T")), 

u^{x,T)ip{x, —,T)dxdT —>■ / / / uo{x,y,T)ip{x,y,T)dydxdT, (6-1) 
Jn e Jo JnJT" 

as e ^ 0. 



We have the following extension of a compactness result first proved by Nguetseng |28] 
and then further developed by Allaire in [29, |3^, 31| and by Holmbom in p^]. 

Theorem 6.1. Suppose that {u^} is a uniformly bounded sequence in L^([0, s]; 

Then there exists a subsequence, still denoted by {u^}, and a function uq = Uo{x,y,T) in 

L^([0, s]; L^(r2 X T")), such that {u^} two-scale converges to uq. 

The relation between u and mq is explained in the next theorem. In fact, by choosing test 
functions which do not depend on y we have the following (see e.g. j29||): 

Theorem 6.2. Suppose that {u^} two-scale converges to uq, where uq = uo{x,y,T) in 
L^([0, s]; L^(r2 X T")), then {u^} converges to uq weakly in L^([0, s]; L^(r2)), where 

u5ix,T) = / uo{x,y,T)dy. 

In other words: 

U = Uq. 



Remark 8. The results of Theorem 6.1 and Theorem 6.2 remain valid for the larger class 
of admissible test functions (p G L'^{[0, s] x $7; C{T"')), see e.g. [^]. 



7 Homogenization 

With the help of the results from Section 4, Section 5 and the two-scale compactness result 
Theorem 6.1 from the previous section we can now state the main results of the paper. 

Theorem 7.1. Consider the Navier-Stokes system (4-2) ■ Suppose that the initial data 
u^{x,Q) = it°(x) and 6^{x,0) = 0^(x) two-scale converge to unique limits u^{x,y) and 
0^{x,y), respectively. Then, as e ^ 0, the following quantities two-scale converge, 



Pe Po, 
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where po = Poi^, t)- Moreover, there exists a function pi = pi(x, y, r) such that 
Vpe VxPO + '^yPl, 

the functions uq, pq, pi and 9q, being the unique solutions to the Navier- Stokes system 
(4-3) with initial data uo{x,y,0) = u^{x,y), 6o{x,y,0) = 9^{x,y) and boundary data T"- 
periodic in the variable y. 

Before we prove the theorem we state a corollary 

Corollary 7.1. Consider the Navier-Stokes system (4.2). The following quantities two- 
scale converge, 

e^^/^tie Uq, 
Pe Po, 

Vpe '^xPO + '^yPl, 
Te ^ To, 

where uq, Po, Pi and Tq are the unique solution to the Navier-Stokes system (4-3). 

Proof of Theorem 7.1. By assumption the initial data admit unique two-scale limits. 

Consider now the first equation of (4.2), 

^ + (ue- V)ue - e^/'^uAue + Vpe = e26e, 'm Q X R+. 
ot 

Let s > and choose test functions (p G Cq°{Q x [to, s]; C°°(r")). By the results of 
Section 5, all the sequences {e^^/^Ug}, {^^}, {Pc} and {9^} are uniformly bounded in 
L^([to, s]; L^(fi)). Therefore, according to Theorem 6.1, they admit two-scales limits. In 
order to identify these limits we multiply each of these terms by the smooth compactly 
supported test function ip{x, f ,t). For the time derivative we get (recall that r = t/-\/e) 

[ [ --^(x,t)ip(x, —,T)dxdT = e~^^'^ [ [ Ue-^dxdr. 
Jo Jn dt € Jo Jq dr 

Sending e ^ and integrating by parts yield, 

ill ^p-^{x,y,r)dydxdT. 
Jo Jn Jt" ciT 

For the second (inertial) term we consider 

/ / [{ue ■ Vue{x, t) - UoVyUo{x, y, T)]ip{x, - ,T)dxdT 
Jo Jn e 

= [ [ [e"^/^'"6 • (e^/^V-Ue - VyUo) + (e-^/^u^ - UQ)VyUQ]ip{x, -,T)dxdT. 
Jo Jn ^ 

By considering VyUoif{x, f ,t) to be a test function in the second term on the right hand 
side this term immediately passes to zero in the two-scale sense. For the first term the 
Schwarz's inequality yields 



/ / e ^^"^Ue ■ {e^^'^Vuc — VyUo) ■ (p{x, —,T)dxdT 
Jo Jn £ 
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< ||e-^/V||||(e^/^Vue - VyUQ)'f{x, ^,r)|| 

<C||(eV2vn,-V,uoMx, V)ll, 

where all the norms are in L2([0,s];L2(f^)). In order to pass to the limit in the right 
hand side we consider the usual L^-moUifications Vu^^^ and VyUQ^^ of e^/^Vn^ and VyUQ, 
respectively. Since the mollified functions pass strongly to e^/^Vug and VyUQ, respectively, 
as /X — 0, we have, for jj, sufficiently small, say /x < /xq and for every y in 

IKe^/^V^Xe - "^yUoM < II (Vm,,^ - VyUoM + S, 

where 6 is arbitrarily small, independently of e. This inequality still holds true if we take 
the supremum in y over T". Thus, for every /x < /xq, the right hand side will tend to zero 
as e tends to 0, by the uniqueness of the two-scale limit V^uo of e^/^Vug. Consequently, 
we have proved that 

-iXeVUg — UoVyUQ 

in the two-scale sense. For the third term we get, by the divergence theorem. 



.^3/2^ 



/ / Au^{x,T)(p{x,—,T)dxdT 
Jo Jn £ 

u / u^{x, T)Ay(p(x, —,T)dxdT + terms tending to zero as e — > 0. 
Jo Jci £ 



Sending e — ^ yields, after applying the divergence theorem again, 

-V I II AyUo(p{x,y,T)dydxdT. 
Jo Jo. Jt^ 

For the right hand side we immediately get 

/ / e29e{x,T)ip(x,-,T)dxdT^ / / / e2do^{x, y,T)dydxdT. 
Jo Jn ^ Jo Jci Jt" 

For the fourth term (the pressure) we have to be a bit more careful. Let us multiply the 
first equation of (4.2) by €ip{x, f ,t). For the pressure term we get 



Jo Jn 



X f x 

Vpe{x,T)ip{x, —,T)dxdT = / / Pf:{x, T){edrvx + divy)ip{x , — , T)dxdT. 
^ Jo Jn e 



A passage to the limit, and an application of the divergence theorem, using the fact that 
all other terms vanish, gives 

/ / / VyPoix,y,T)(f{x,y,T)dydxdT = 0, 
Jo Jn Jt" 

which implies that po does not depend on y. We now add the local incompressibility 
assumption on the test functions ip, i.e. divy(p = and multiply the pressure term by 
ip{x, f ,t) and apply the divergence theorem, 

/ / Vp^{x,T)(p{x, —,T)dxdT = / / p^{x,T)diYx(p{x,—,T)dxdT. 
Jo Jn £ Jo Jn ^ 
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A passage to the hmit, and an apphcation of the divergence theorem, gives 

/ / / VxPoix,T)(p{x,y,T)dydxdT. 
JO JnJT" 

Collecting all two-scale limits on the right hand side gives 

/ if ~ 'a (■"o ■ '^y)'^o + ^\uo - VxPoj^dydxdr = 0. 

Jo Jn Jt" (^'^ 

Since divycp = we can argue as in Remark 2 and conclude that there exists a local 
pressure gradient VyPi(x, y, r) given by 

^yPiix, y,T) = f - - (no • Vj,)no + vlS^yUQ - VxPo- 

Let us now consider the second equation of (4.2). We already know that the sequence {6^} 
is uniformly bounded in i^^([0, s\;L?{p)). We multiply by e^/'^ip as above and for the time 
derivative we get 



By letting e ^ we get 



6f:-^dxdT — / / / —^ipdydxdr. 



/o Jn 9t Jq Jq dr 

For the non-linear term we have 



eV2 



/ {ue ■ V)0^ipdxdT 
Jo Jci 

I I {Ue-Vx)<^OedxdT + t'^/'^ [ [ {u^ ■ V y)ip6edxdT. 

Jo Jn Jo Jn 



The first term on the right hand side immediately passes to zero. For the second term we 
argue as above and consider the difference 



/ / ((e-^/^Ue ■ Vy)ip9e - (no • Vy)ip9o)dxdT 
Jo Jn 

= 11 ((e"^^^^i6 • yy)^Oe - {UO ■ Vy)(p9e)dxdT 

Jo Jn 

+ / / ((^*0 • ^y)'pGe - {uo ■ Vy)(p9Q)dxdT. 

Jo Jn 



By considering {uq ■ Vy)^ to be a test function in the second term, this term immediately 
passes to zero in the two-scale sense. For the first term we get, by the Schwarz's inequality, 

/ / {{e-^'^u, ■ Vy)ip9e - (no • Vy)ip9e)dxdT < C\\{e-'/\e - uo) ■ Vyifl 
Jo Jn 

where the norm is in L^([0, s]; L^(r2)). We introduce, as above, mollifiers and consider 
the sequence {n^^} which converges to e~^/^ng strongly as — > 0. Arguing as above, we 
choose sufficiently small to get 

||(e~^/^ne - no) • Vyip\\ < [Kn^,^ - no) • Vy(p\\ + 6, 
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where 5 is arbitrarily small. Thus, by sending e — > 0, 

in the two-scale sense. For the third term we get, by the divergence theorem, 

—e^K I / A9f(pdxdT 
Jo Jci 

= —K / 9eAy(pdxdT + terms tending to zero as e — 0. 

JO Jn 

We let e ^ and get, by the divergence theorem, 



Jo Jci 



yipdxdr — > —K III AyOofdydxdT. 
JciJt" 



The right hand side of the second equation in (4.2) will vanish since is of order e^/^. By 
Theorem 4.2 the system (4.3) has a unique solution {uo,po,pi,9o} and, thus, by unique- 
ness, (4.3) is two scales homogenized limit of the system (2.1). Also, by uniqueness, the 
whole sequence converges to its two-scale limit and the theorem is proven. ■ 

8 The mean velocity field 

In this section we derive the mean field uo for the velocity. Let us consider the first 
equation of (4.3) 

+ {uo ■ Vy)uo - I'AyUo + VyPl = CnOo - ^xPO, 

divyUo = 0, divx( / m dy) = 0, 

in X T" X R+ with r"-periodicity as boundary data in y. 

By letting K = K{y,y',T) denote the heat kernel we can write the solution to the first 
equation of (4.3) as an integral, 

uo{x,y,T) = / K{y,y',T)uoix,y',0)dy' 

+ / / K{y,y',T- s){en0o{x,y',s)-VxPo{x,s)-VyPi{x,y',s) 
Jo Jt^ 

-{uq ■ Vy)uoix, y', s))dy'ds. 

Now, by the decay of the heat kernel, for r sufficiently large, the first term becomes 
arbitrarily small. An averaging of the second term over in y gives 



uq{x, t) = I {cnOo - VxPo - (uo ■ Vy)uo){x, s)ds. 
Jo 



10 

The divergence theorem gives 



/ {uq ■ Vy)uo dy = - {diYyUo)uody + [uq ■ %)mo dSy, 

JT" Jt" JdT^ 
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where Uy is the local unit normal and Sy the local surface element. Now, by the local 
incompressibility we get 

{divyUo)uody = 0. 

'Jpn 

Moreover, for the boundary integral we get, in the two-dimensional case, 

/ UiniUodSy= -U2Uo{yi,0)dyi+uiUo{l,y2)dy2 
JdT2 Jo 

+ / U2Uo{yi,l)dyi - uiuo{0,y2)dy2 = 0, 
Jo 

by the T^-periodicity of uq = {ui,U2)- Consequently 



{uo ■ Vy)uo = 0. 

The computation in the three-dimensional case is similar. Therefore the mean field reduces 
to 

u^{x,t)= / {enOo -VxPo){x,s)ds. (8.1) 
Jo 

We can now use the incompressibility, which when applied to (8.1) gives 

dw^{enOoix, t) - V:r.po{x, t)) = 0. (8.2) 
This says that 

e^o{x, t) - VxPoix, t) = H{x, r) (8.3) 

where if is a divergence free (rotational) field. 

The field H can be determined explicitly by solving the global equation (8.2), (in x), for 
the pressure po- For this purpose one needs to impose the appropriate boundary condition 
on the pressure in order to close the system. We impose the Neumann condition 

n ■ VxPo = 

for the pressure. In fact Lions [45| and Sanchez-Palencia [^] impose the condition 

n • Mo = 
which when inserted in (8.1) gives 

• V^po = 0. 
We consider 

AxPo{-,t) = dw^enOoi-jT), in Q, 

(8.4) 

n ■ VxPo{-,t) = 0, on dQ, 
for every r > 0, and obtain 

H = 7T{eJo) = -V X (A-i(V X eJo)), 
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where 7r(en^o) denotes projection onto the divergence free part of the conduction term. 
Collecting the results from the above discussion we can now express the mean field as 

t f^/^ 
uo{x,—p)= / TT{en0o){x,s)ds. 
Jo 

This gives the contribution of the conduction and the global pressure to the small scale 
flow. We have let the local flow settle down and averaged (in y), denoted by overbar, over 
the unit ceU T", n = 2, 3. 

Remark 9. If we take boundary layer effects into account, we can, considering the sim- 
plest case of a boundary layer, specify the value of the global pressure gradient VxPo at 
the boundary, of the boundary layer. For instance we can put 

n ■ V^po = f. 

This will result in the additional term G * / in the mean field, where G is the usual 
Neumann kernel for the Laplacian. In this case the mean field becomes 

uo{x,—f)= / {■^{enOo) +G * f) {x,s)ds. 



Remark 10. The above formulas give the mean field fiow with the global convection 
(rolls) taken out. This was done by imposing the periodic boundary conditions on the 
local cell, see the introduction. The two-scale convergence carries over to a larger class 
of test functions which are non-periodic , see [48|. This will not be repeated here, but 
to be able to compare the true (collective) mean field we compute the average over the 
convection term in the non-periodic case. By the Taylor and mean value theorems we get, 
in the two-dimensional case, 

/ UiniUQdSy= {-U2Uo{yi + Yi,Y2) + U2Uo{yi+Yi,l + Y2))dyi 
JdT'^ Jo 

+ / {uiuoi'i- +Yi,y2 + Y2) - uiUo{Yi,y2 + Y2))dy2 
Jo 

= U2 [ {uo{yi + Yi,Y2) - uo{yi + Yi,l + Y2))dyi 
Jo 

+ui [ {uo{l+Yi,y2 + Y2)-uo{Yi,y2 + Y2))dy2 + 0{A) 
Jo 

dun dun ^ , » , „ ^ / ^ n 

= U2^ + ui-^ + 0{A) = uo-Vu^ + OiA), 

Ol2 Oil 

where A = {yi — yi,y2 — ^2)) ^i^y2) is the point in where u attains its mean value, 
and (11,12) is the location of the box. The error term 0(A) < (7||uo|||A| is bounded by 
Theorem 4.2. In global coordinates we therefore get 



no • Vuo = Mo • Vuo + 0(e), 
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and similarly 

If we insert this into the expression for the mean field and differentiate with respect to the 
fast time variable r, we obtain, as e — > 0, 

duQ _ ^ , , 

— + uo • Vuo + Vpi = 7r(e26'o), 

i.e. a forced Euler equation, where uq is a function of the scaled variables (r, Y) = 
(t/-y/e, x/e). The computation in the three-dimensional case is similar. Thus the local flow 
satisfies the Navier-Stokes equation, whereas the mean flow satisfies the Euler equation. 
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